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Graph-Directed Sprays and Their Tube Volumes via 
Functional Equations 

Derya QELiK * §ahin KOQAK * t Yunus OZDEMiR* 

A. Ersin UREYEN * 


Abstract 

The notion of sprays introduced by Lapidus and his co-workers has proved useful 
in the context of fractal tube formulas. In the present note, we propose a more 
general concept of sprays, where we allow several generators to make them more 
convenient for applications to graph-directed fractals. Using a simple functional 
equation satisfied by the volume of the inner e-neighborhood of such a generalized 
spray, we establish a tube formula for them. 
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1 Introduction 

Tube volumes for fractals has been an interesting topic of research in the last decade. 

M. Lapidus and his coworkers established several tube formulas as series in terms of 

residues of some associated zeta functions |7]. Recently, we have proposed a simple 

alternative approach to tube formulas for self-similar sprays via functional equations 

([!]). The so-called sprays introduced by Lapidus and Pomerance in [10] (to be defined 

below) are a convenient notion in dealing with fractal tube volumes. 
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Beyond the self-similar fractals, the next family of interest is the class of graph- 
directed fractals. For tube volumes of graph-directed fractals there are formulas (see 
m, which are analogous to the tube formulas of Lapidus and Pearse ([8]). So far as 
we know, the notion of spray has not been extended to the graph-directed setting. In 
the present note, we want to propose the notion of graph-directed sprays and establish 
a tube formula for them via functional equations. 

We now first recall the notion of spray in Euclidean space. Let G C M n be a bounded 
open set. A spray generated by G is a collection S = (Gi)^ of pairwise disjoint open 
sets Gi C M n such that G* is a scaled copy of G by some Aj > 0; in other words, Gi is 
congruent (i.e. isometric) to A iG. The sequence (Aj)jgN is called the associated scaling 
sequence of the spray. In applications one has typically A* < 1 and often Ao = 1 so that 
Go is equal (or isometric) to G. Furthermore, it is meaningful to assume A” < oo to 

iSN 

make the volume of |^J G* finite, as one deals with the inner tube of |^J Gj (the inner 

ieN ieN 

e-tube of an open set A C M n is the set of points of A with distance less than e to the 

boundary of A). 

The most important class of sprays is the self-similar sprays, for which the scaling 
sequence is of a very special type. 

Let {?’i, r 2 ,..., rj} be a so-called ratio list (i.e. 0 < rj < 1 for j = 1, 2,..., J) and 
consider the formal expression 
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If the scaling sequence (A;)^ °f a spray is given by the terms of the series on the 
right-hand side of © for an appropriate ratio list, then the spray is called a self-similar 
spray. 

To give a flavour of tube formulas, we note the following theorem (0, [9j, 0). 


Theorem 1 Let (Gj)jgN be a self-similar spray generated by G C M n with a scaling 
sequence associated with a ratio list {ri,r 2 , ■ ■ ■, rj}. Assume G C M n to be monophase, 
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i.e. let its inner e-tube volume function Vg(b) be given by 


V G (e) 


' n— 1 

]T K t S n - 

< i=0 

Vol(G) 


0 < e < g 

£>g, 


(2) 


where g is the inradius of G (i.e. the supremum of the radii of the balls contained in 
G). Then the volume V\j ^(e) of the inner e-tube ofUGi is given by the formula 


VuGi(z) = res{C{s)e n s -u) , for e < g , 

wG2)U{0,l,2,...,n-l} 


n 

where C(s) = 1 _ (rf+ 1 .„ +r , ) 

i =0 

roots of the equation rf + • • • + r j 


s — i 


: Hi 


(with K r 


= 1. 


—Vol(G)) and D is the set of complex 


The motivation behind the notion of the self-similar spray is that they naturally 
emerge as the “hollow spaces” in self-similar fractals. For example, if you start with an 
interval, and construct a Cantor set by deleting successively the open middle thirds of 
the intervals, the collection of deleted open intervals constitute a 1—dimensional self¬ 
similar spray with a scaling sequence associated with the ratio list the generator 

being the first deleted middle third. Likewise, if you start with a triangle and succes¬ 
sively delete the open middle fourths of the triangles to obtain in the end a Sierpinski 
Gasket, then the collection of the deleted open triangles constitute a 2—dimensional 
self-similar spray with a scaling sequence associated with the ratio list {|, ^}, the 

generator being again the first deleted middle fourth. 

One can in principle allow any scaling sequence for a spray, but to obtain manageable 
tube formulas some sort of restrictions seem (as to yet) to be necessary. To be associated 
with a ratio list is, for example, such a condition. A weaker condition (called “subshift 
of finite-type”) was formulated to handle the hollow spaces of graph-directed fractals 
in [3], but a more natural approach demands a consideration of sprays with more than 
one generator, since the hollow spaces of graph-directed fractals are composed of copies 
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of several generators (see Figures I3ll5lb 

In Section 2 we propose this more general concept of sprays with several generators. 
In Section 3 we consider a natural functional equation for inner tube volumes of graph- 
directed sprays, formulate a multi-dimensional renewal lemma to handle it and establish 
an inner tube formula for graph-directed sprays as our main result (Theorem [2]). In 
Section 4 we give the proof, taking into account the additional difficulties arising from 
the presence of the Mauldin-Williams matrix. We give also an explicit sufficient region 
for e on which the main theorem holds. 

2 Graph-Directed Sprays 

Let Q = (V,E,r) be a weighted directed graph with weights r : E —>• (0,1). For an 
edge e € FI, we denote the initial vertex of e by i(e) and the terminal vertex by t(e). 
For vertices u,v € V, we denote the set of edges from u to v by E uv and the set 
of edges starting from u by E u . If E u 0 for all u € V, such a graph is called a 
Mauldin-Williams graph. If any two vertices u and v can be joined by a (directed) 
path, then the graph is said to be strongly connected. We will generally assume that 
the Mauldin-Williams graphs be strongly connected. 

We define the weight of a path a = eie 2 • • • e*, by r(a) = r(ei) • r(e 2 ) ■... • r(efc). a is 
called a path from the vertex u to v if z(ei) = u and t{ejf) = v. We also write i(a) = u 
and t(a) = v. We assign an empty path 4> u to every vertex u with weight r{(f u ) = 1. 
Now we define graph-directed sprays. 

Definition 1 Let Q = (V, E, r) be a Mauldin- Williams graph and G u ( u E V) be 
bounded open sets in R n . A graph-directed spray S associated with Q and generated 
by the open sets G u (u £ V) is a collection of pairwise disjoint open sets G a in M n 
(where a is a path in the graph), such that G a is a scaled isometric copy of G t ^ with 
scaling ratio r(a). 


Remark 1 Note that G u , with u € V, is a generator and G a , where a denotes a path, 
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is a copy of a generator. If a = cj> u then G^ is a scaled isometric copy of G t u n ^ = G u 
with scaling ratio r(4> u ) = 1, i.e. G^ is an isometric copy of G u . 

Remark 2 Note that if Q has only one node this notion reduces to the ordinary notion 
of a self-similar spray generated by a single open set with a scaling sequence associated 
with the ratio list consisting of the weights of the loops around the single vertex. 

The spray S can naturally be decomposed into subcollections 

Su = { G a | i(a) = u }. 

Notice that the subcollection S u is also composed of pairwise disjoint scaled copies of 
all generators G v with scaling ratios r(a) for paths a starting at u. This decomposition 
of S into the subcollections S u will prove useful in establishing tube formulas. 

The motivation for the definition of graph-directed sprays comes, in analogy to the 
motivation of self-similar sprays, from hollow spaces of graph-directed fractals. Let us 
very briefly recall the notion of graph-directed fractals. 

Let Q = (V, E, r) be a Mauldin-Williams graph, (A u ) u e y be a list of complete 
subsets of M n and let f e : A t ( e ) —>• A ^ be similarities with similarity ratios r(e). Such 
an assignment is called a realization of the graph Q in R n . Given such a realization, 
there is a unique list ( I\ u ) u& y of nonempty compact sets with K u C A u (u € V ) 
satisfying 

R u= U U /e(^) 

v&V e&E 

UV 

for all u £ V (|5])- 

In favorable cases the maps f e can be restricted to the convex hull of the graph- 
directed “attractors” K t / e ^ and these attractors can be imagined to be formed by delet¬ 
ing successively pieces of the convex hull analogous to the construction of the Cantor 
set or the Sierpinski Gasket by deleting successively pieces of an interval or a triangle. 
The collection of deleted open pieces will constitute a graph-directed spray in the above 


Graph-Directed Sprays and Their Tube Volumes via Functional Equations 


6 


defined sense, with so many generators as there are nodes of the graph. Before making 
this idea precise, it will be best to study an example. 

Example 1 Consider the Mauldin-Williams graph with V = {1,2}, with 9 edges and 
the corresponding weights as shown in Figure [/} 



Figure 1: A Mauldin-Williams graph with 2 nodes and 9 edges (the weights are shown 
in parenthesis). 

Let A\ and A% be the square and the triangle in M 2 as shown in Figured (a)- The 

similarities associated with the edges are indicated in Figure d (b). The graph-directed 

fractals K\ and K 2 of the system are shown in Figured (A\ \ K\) U ( A 2 \ K 2 ) is a 

collection of connected open sets which constitute a graph-directed spray with generators 

G\ and G 2 satisfying Definition [7j (see FiguresThe generators G\ and G 2 thereby 

are defined by G u = A° u \ / e (A t(e) ) for u = 1,2. 

e€E u 

Note that the subcollections 5i = {G a \ i(a) = 1} and S 2 = {G a \ i(a ) = 2} contain 
copies of both of G\ and G 2 so that we can not view either of them as a spray in the 
classical sense. 

To illustrate the formation of G a for a path a, we give several examples in Figures 0 
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Figure 2: A realization of the Mauldin-Williams graph of Example [Tj (shown in 

Figure [T]). (a): The complete spaces associated with the 2 nodes, (b): The similarities 
associated with the edges. 
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Figure 3: The attractors of the realization shown in Figure [2] of the Mauldin-Williams 
graph of Example [l] 


( 0 , 1 ) 


( 2 , 1 ) 


( 1 , 0 ) 
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Figure 4: The generators of the realization shown in Figure[2]of the Mauldin-Williams 
graph of Example [TJ 




Figure 5: The hollow spaces of the realization shown in Figure [2] of the Mauldin- 
Williams graph of Example [□ (<Si left, £2 right). 
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path a 


a 


a 



Figure 6: A few examples of actions of paths on generators for the realization shown 
in Figure [2] of the Mauldin-Williams graph of Example [TJ 
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From the point of view of our present concern to establish inner tube formulas for 
graph-directed sprays, the special positions of the scaled copies of the generators are 
not important as long as they are pairwise disjoint. But if one wishes to compute tube 
volumes of fractals, one should be careful in relating the tube of the fractal to the inner 
tube of an associated spray. In the above example the e-tube volume of K\ can be 
expressed as the sum of the inner e-tube volume of the collection <Si and the outer 
e-tube of the square A\ (likewise for K'f) (see Figure 0. But this happy relationship 
does not hold always as the following example shows. 




Figure 7: The e-tubes of the attractors of Example [T| as union of the inner e—tube of 
the hollow spaces and the outer e—tube of the convex hull of the attractor. 


Example 2 As the Mauldin-Williams graph we choose the same graph in Figure Q] 
with the only difference that we delete the edge e\ 2 - We choose the same realization of 
this graph in M 2 discarding the map corresponding to e\ 2 . The emerging graph-directed 
fractals L\ and L 2 and the corresponding graph-directed spray with generators H 1 and 
H -2 satisfying Definition [7] are shown in FigureslMB 




Figure 8: Attractors of Example [2] 
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Figure 9: Generators and the hollow spaces of Example [2] showing that the e— tube 
of the attractors needn’t be the union of the inner e— tube of the hollow space and the 
outer e—tube of the convex hull of the attractor. 

In this example the e-tube volume of the graph-directed fractals can not be mean¬ 
ingfully related to the inner e-tube volume of the graph-directed spray. Two types of 
problematic boundaries of generator copies are indicated in Figure [21 by dotted lines. 

The simple relationship between the fractal tube volume and the inner spray volume 
observed in Example [T] still remains true for a more general class of graph-directed 
systems, if the following assumptions hold: 

i) dim(C' u ) = n, where C u is the convex hull [ K u \ of K u . 

ii) (Tileset condition) The open set condition should be satisfied with O u = C°. We 
recall that (see 0 ) the graph-directed system satisfies the open set condition if 
there exists a list ( O u ) u& v of open sets O u C M n such that, 

(a) for any e € E uv , f e (O v ) C O u , 

(b) for any two distinct ei,e 2 € E u , f ei (O t (e i)) Fl /e 2 (Ot( e 2 )) 

iii) (Nontriviality condition) C° ^ U /.(<?«.))• 

e£E u 
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iv) (Pearse-Winter condition, [12]) 8C U C K u . 

Now, if we define G u = C° \ |^J / e (G(e)) then we get a graph-directed spray S 

e(E.E u 

generated by the open sets ( G u ) ueV with G a = G eie2 ... ek = f ei fe 2 ■ ■ ■ fe k (G t ( efc )) for a 
path a in the graph Q. 

Under the above conditions one can compute the e-tube volume of the graph- 
directed fractals with the help of the inner tube volume of the graph-directed spray 
as in Example |T] So we now consider the inner tube volumes for graph-directed sprays 
in the next section. 

3 Inner Tube Volumes of Graph-Directed Sprays via Func¬ 
tional Equations 

Let Q = (V, E, r) be a Mauldin-Williams graph, ( G u ) u& v bounded open sets in M n and 
S be a graph-directed spray associated with Q and generated by the open sets (G u ) ug y. 
Let S u be the subcollection of the spray S corresponding to the paths with initial vertex 
u£V. 

The volume of the inner e-tube of the collection S u satisfies the following functional 
equation (for all u E V) as one can easily verify: 

v&.(£) TE'J %„( £ /r,) + VaM. (3) 

veVeeE 

uv 

Our strategy will be, as in the self-similar case (0), to apply the Mellin transform 
to this functional equation and then try to recover the volume function by applying the 
inverse Mellin transform. To apply the Mellin transform we need an estimate of Ys u (s) 
as e —> 0. We now formulate a multi-dimensional renewal lemma which will enable us 
to find such an estimate (For the one-dimensional renewal lemma see Ti i). 

We recall that for a strongly connected Mauldin-Williams graph the spectral radius 
of the matrix ^4(s) = [a u „(s)] M V&v with a uv (s ) = rf takes the value 1 for a unique 

edzEuv 
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so > 0, which is called the sim-value of the graph and which we denote by D below 
(see HD- We will always assume D < n. 

Remark 3 The assumption D < n is in fact equivalent to the condition that the total 
volume of the graph-directed spray is finite. To see this, one can easily verify that the 
volumes of the subcollections S u can be expressed as follows: 

[Vol(S u )\ U £v = {I + A(n) + A 2 (n) + • • •) [V ol(G u )] u ^y 5 

where [Vol(S u )\ ue y is a column vector. Note that, the matrix power A k (n ) codes the 
contribution of paths of length k to the total spray volume. 

One can now prove the above claim taking into account that the spectral radius 
of t4(s) is a strictly decreasing function for s > 0 Jfff, and for any matrix A, each 
entry of A k tends to 0 as k —>• oo if and only if the spectral radius is less than 1 | 6 j 
Theorem 5.6.12]. 

Lemma 1 Let Q = ( V,E,r ) be a strongly connected Mauldin- Williams graph and 

hu{x) = ^2 ^2 r e h v(x ~log(l/r e )) + fi u (x) (ueV), (4) 

vGVeGE 

uv 

be a system of renewal equations on R, where D is the sim-value of the graph. Assume 
'fiu(x) = 0(e~ T \ x \) for some r > 0. Let {h u ) u( zy be a solution of this system of renewal 
equations. If h u tends to 0 for x —» —oo for all u € V, then h u is bounded (for all 
u € V). 

Proof. Let 7 = min 

e&E 

Since the Mauldin-Williams graph Q is strongly-connected, the corresponding Mauldin- 
Williams matrix ^4(s) is irreducible and by the Perron-Frobenius theorem, for s = D 
the spectral radius 1 is also an eigenvalue with a positive eigenvector p = (p u )u&v (with 
p u > 0 ) so that we have 
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Pu = J2 a uv (D)p v = ^2 r e Pv ( 5 ) 

veV veV e£E uv 

Since the functions h u ( u E V) tend to zero for x —» —oo, one can choose xq E M 
such that \h u (x)\ < p u for x E (— 00 , Xo]. 

Let x E [xq , xq + 7 ]. From Q and (|5]1 


\h u (x)\ < E E rf p v + sup \ip u (x)\ 

veVeeE uv xe[a:o,0:0+7] 

< p u + sup I tpu(x)\. 

x£[xo,xo+"f\ 

By the assumption on ij) u , we can find an M such that \ip u (x)\ < p u M e ~ T ^ (for 
all u G V). Hence for x G [xo, xo + 7 ] 

\K(x)\ < p u ( 1 + M sup e _r l 3: l J . (6) 

\ o:G[o; 0 , 0:0+7] J 

Since \h u (x)\ < p u for x E (— 00 , Xo], the inequality (JUJ) holds for all x E (— 00 , Xq + 7 ]. 
Now, let x E [xq + 7 , xq + 27 ]. As above, 


I+MI fi E E rj?p v ( 1 + M sup e 


-t\x\ 


v£V e£E u 


x£[xq ,XQ-\-'y] 


+ SUp \lp u (x)\ 

xe [ 0 : 0 + 7 , 0 : 0 + 27 ] 


< p u ( 1 + M sup e T \ x \ ] +p u M sup e 


-t\x\ 


X£ [0:0,0:0+7] 


= Pu 


xe [0:0+7,0:0+27] 
-t\x\ 


1 + M sup e~ T ^ + M sup e 

x€[x 0 , 0:0+7] o:G [0:0+7,0:0+27] 


The above inequality clearly holds for all x E (— 00 , xo + 27 ]. Repeating the above 
argument, we see that for all x E M, 


\K(x)\ < p u i 1 + M ^2 sup 


T |as| 


k=0 o;G[x 0 +fc7,o:o+(A:+l)7] 


This shows that h u is bounded on M. ■ 

Now we can derive an estimate for Vs u (e) as e -+ 0. 
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Lemma 2 Assume, there exists an a > 0 such that Vg u (e) = 0(e a ) for e —» 0. We 
further assume that n — a < D < n, where D is the sim-value of the graph. Then, it 
holds Vs u (e) = 0(e n ~ D ) as e —>• 0. 

Remark 4 If the generator G u is monophase or pluriphase (i.e. the volume of the 
inner e-tube of G u is piecewise polynomial), then a > 1 and the assumption n — a < D 
reduces to n — 1 < D. In applications to graph-directed fractals, G u will be mostly 
monophase or pluriphase, but in graph-directed sprays G u can be themselves of a fractal 
nature, so that a could be less than 1 . 

Vc (V) 

Proof of Lemma O Let us define Ws (e) = —ft to obtain 

/ £ n —D 

Ws* (e) = E E r ? Ws ” (e/re) + 

v(z.V e(E.E uv 

from the equation (f3|) . Let us now apply the change of variable e = e~ x in the equation 
(HD- We obtain the following system of renewal equations on M 


h u (x) = Y E r e hv{x - logl/r e ) + ipu(x), 


veV e£E u 


where h u (x) = Ws u (e x ) and tf u (x) = 


V Gu (e~ x ) 


g—x(n—D ) 

Now we have to verify the assumptions of Lemma [lj 

V s (e) 

hu(x) = h u (— loge) = tends to zero as x —>• —oo, i.e. e — > oo, since D < n 

and the volume of the spray is finite. 

To check the assumption on if u , choose r = min{n — D, D — n + a}. 

By Lemma [TJ h u {x ) = Su ^J is bounded, so that h 5 u (e) = 0(e n ~ D ). ■ 

e n u 

We will now apply the Mellin transform to the equation (j3]) and to this end, it 

V s (e) 

will be convenient to define the auxiliary functions / u (e) = —^—, for u € V (these 
functions can be viewed as a kind of “normed” tube volumes). The system ([3]) of 
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functional equations transforms into the following system: 

/«(£) = E E r-• (8) 

v€V e&E uv 

Recall that the Mellin transform of a function / : (0, oo) —> M is given by 


~ r°° 

M(f)(s) = f(s)= f(x)x s ~ 1 dx. 

Jo 


If this integral exists for some c € M and if the function / is continuous at x E (0, oo) 
and of bounded variation in a neighborhood of x, then f(x) can be recovered by the 
inverse Mellin transform (I f3|) 


- lim 

27ri T—>oo 


rc+iT 

Jc-iT 


/(«) 


x S ds. 


The function f u is continuous and f u (s) = 0(e ~ n ) as e —>• oo. If for some a > 0, 
Vq v (e) = 0(e a ) as e —>• 0 and n — a < D < n, then by Lemma O f u (e) = 0(s~ D ) as 
e —>• 0. So the integral 

POO 

/ /«(e)e s_1 de' 

Jo 

exists for any s with D < Re(s) < n. Likewise, the integral 


f 

Jo 


e 3 - 1 de 


exists for n — a < Re(s) < n. We can then take the Mellin transform of (J8]) to obtain 

w*) = EfE <)/»('» 

vev \eeEuv / v 7 

for D < Re(s) < n. This system can also be written as a matrix equation 


F{s) = A(s)F(s) + $(s) (D < Re(s) < n ), 
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where F(s) is the column vector 
is the Mauldin-Williams matrix. 



ugV 1 


$(s) 



(«) 


u&V 


and 7L(s) 


Lemma 3 For Re(s) > D, the matrix I — A(s) is invertible, so that it holds 


F(s) = [I-A(sT 1 <h(s) 


(9) 


for D < Re(s) < n. 

Proof. This is a consequence of some well-known results from matrix algebra. For s € 
M, s > D, the spectral radius p(A(s)) is less than 1 ([5j). Then by [6], Theorem 5.6.12] 
lim A k (s) = 0 entry-wise. For arbitrary s € C with Re(s) > D, we have 

k —^OO 


.mi = 


£ 

G^Euv 


< 


rf e{s) = a uv (Re(s)). 


e(HEu 


This holds for the entries of A fc (s) and A fc (Re(s)) also, giving lim A k (s) = 0 entry- 

k —^oo 

wise. We then have by [6] Theorem 5.6.12], p(A(s)) < 1, and thus I — A(s) is invertible. 


We can write the matrix equation |9]) also as follows: 



uev 


1 

det(I — A(s)) 


[adj(7 — -A(s))] w 



w 


uev 


or 


AM = dot(/ * A(s)) E M (VdAj („) for all UEV. (10) 

We now apply the inverse Mellin transform to the equation (1101) : For D < c < n, 


fu(e) 


- f 

2vri J c _ 


— y 

2 iri ' 


e+ioo 


fu{s)e S ds 


vev 


pc+ioo 


adj(J - A(s)) w 
det(J — .A(s)) 


■ M 



(s)s s ds 
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Definition 2 Let Q = ( V,E,r) be a Mauldin-Williams graph, G u (u € V) bounded 

open sets in and S a graph-directed spray associated with Q and generated by the 

open sets G u . Let A(s) = veV with a uv (s ) = r s e be the Mauldin-Williams 

e(zE uv 

matrix of the graph. We define the geometric zeta function of the graph-directed spray 
with respect to the node u E V as follows: 


c«w = X M () w , 


veV 


det (I — A(s)) 


for D < Re(s) < n where D is the sim-value of the Mauldin-Williams graph. 


/ u (e) can now be expressed as 


/«(e) 


1 

27ri 



e S ds. 


At this point, we need some assumptions about the inner tube volumes of the 
generators to manipulate this expression further. We assume that the generators are 
monophase or pluriphase. Note that in this case the geometric zeta function which 
is analytic for D < Re(s) < n can be extended meromorphically to the whole plane 
C, as can be seen from the explicit expressions for the Mellin transforms given in the 
Remarks [MSI below. 

We can now express our main result as follows: 


Theorem 2 Let Q = (V,E,r) be a Mauldin-Williams graph, G u [u € V) bounded open 
sets in M n and S a graph-directed spray associated with Q and generated by the open 
sets G u . Let A(s) be the matrix of the Mauldin-Williams graph. 

We assume the generators G u to be monophase or pluriphase. We furthermore 
assume that the sim-value D of the Mauldin-Williams graph satisfies n — 1 < D < n. 

Then for small e, the volume of the inner e—tube of the graph-directed spray S can 




Graph-Directed Sprays and Their Tube Volumes via Functional Equations 


19 


be expressed pointwise as the following residue formula: 


V s {e) = J2 X res(Cn(s)e n s ;w), 

ueV we®U{0,l,2,...,n-l} 

where D is the set of zeros of det(I — A(s)), which we call the complex dimensions of 
the graph-directed spray (For an exact bound for e, see the last paragraph of the proof). 

Remark 5 If the generator G u is monophase with tube formula 


( n— 1 


Vcjs) = { 


Y,< £n ~ l > 0 <e<g u 

i= 0 

Vol(G u ) , £ g u , 


then 


where = —Vol(G u ). 


M 


V Gu {e) 


w = X K 

i=0 


s — x 


Remark 6 If the generator G u is pluriphase, let us assume that it has the tube formida 


V Gu (s) = { 


X k T ,U £n 1 > 9rn-hu < £ < 9rn,U , m = 1 , 2 ,..., M u 

i =0 

Vol(G u ) , £>9u, 


where go,« = 0, 9m u ,u = 9u (9u the inradius of G u ) and kI ’ 11 = 0. It will be more 
convenient to write above formida as 


Vg u (£■) = X k T' U £n 1 f° r 9rn-i,u <£< 9m,U , rn = 1, 2,..., M u + 1, 

%=0 

where we set = 0 /or i = 0, 1 ,... , n— 1 , k^ u+1,u = Vol(G u ) and 9 m u +i,u = oo. 

Then 

m (w = x x ( / C ,U - «r +1,w> ' 


ra=l 


S — l 
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Example 3 (Example [l] continued) The volumes of the inner e —neighborhood of 
G\ and G 2 are given by the following functions: 


V Gl (e) 


4\/2e - 4e 2 

< 

2 




{ 2± fi) £ ~ (3 + 2y/2)e 2 , £ 


Vc 2 (e) = < 

The corresponding Mauldin- Williams matrix of the graph is 


< 

£ > 


2-V2 

4 

2-^2 

1 


Ms) 


0 

1 

Y s 




and the sim-value of the graph is D = log 2 • The complex dimensions of the 

graph-directed spray are given by 


log 2 


V29 +1 


+ ikp | k G Z > U < log 2 


^ 1 ) +i (k + \ )p|fcez 


27T 

where p = -. Using Theorem [H we obtain the volume of the e—neighborhood of S as 

In 2 


4 28 2 3 

Vs(£) = — £ 2 -— \/2e + Si + S 2 + —(3 + 2\/2)£ 2 — —(2 + \/2)e + S 3 + E4 

7 5 7 5 


where 


Si = 


g.2—(D+ifcp) 

\/291n2 


2\/29 — 2 


7 


-4 


/ y^2 \_DH-ifcz? ( V^ \D—l+iky 

1 2 ’ + 4\/2 > 2 J 


D + ikp 


- 2 


( vd^j-D—2+ifcp' 


D — 1 + i&p D — 2 + ifcp 


+4 -(3 + 2\/2) 


^ 2—\/2 ^D-)-ifcp ^ 2—\/2 )_D—l+ifcp ^ ^ 2— V2 ^D— 2+ifcp ' 


-D + ifcp 


+ 


D — 1 + i/cp 


D — 2 + ikp 
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So = 


r 2-(D'+i(fc+i)p) 

\/29 In 2 
2\/29 + 2 


7 


-4 


^\/2^D'+i(fe+i)p 

-O' + i(/c + \)p 


+ 4^2 


' \/2 \D'— l+iffc+llp 


- 2 


' V2^D'- 2+i(fc+l)p 


D' — 1 + i(fc + D' — 2 + i(fc + 


-4 —(3 + 2\/2) 


^ 2—y/2 '^D r -\-i{k-\-\)y 2 _|_ ^2 ( 2—\/2 ^D'— | ^ 2—\/2 ^D' — 


D' + i (k + \)p 


+ 


D' — 1 + i (k + \)p 8 D' — 2 + i(/c + \)p 


S 3 = 


e 


2 —(D+ikp) 


\/291n2 

v^2 ^-D+ifcp 


-4- 


+ 



z'VIl+ifcp t V 2 \D— 2 +ikp \ 

+ 4v /2^J-— -2^-— 

D — 1 + i kp D — 2 + i/cp / 


^ 2— x/2 ^D+ifcp 2 | v ^2 ( V2 )_D—1+ifcp ^ 2— \/2 ^D—2+ifcp^ ' 


D + ifcp 


D — 1 + i&p 8 D — 2 + i&p 


S 4 = — 


r 2-(D'+i(fc+i)p) 


-4 


\/291n2 

^)L»'+i(fc+|)p 

D' + i(fe + i)p 


(y/ 2 \D'~ l+iffc+4)p 

+ 4\/2 2 


- 2 


( v / 2^D / -2+i(fc+|)p 


V29 + 1 


-(3 + 2\/2) 


D' — 1 + i(fc + |)p D' - 2 + i(fc + y 

^ 2—\/2 2 _(_ y /2 ^ 2—\/2 ^ D '—^ ^ 2 — y /2 ^ D '— 


D' + i (k + h)p 


+ 


D' — 1 + i(/c + 8 D' — 2 + i(fc + i)p 


4 Proof of the Main Theorem (Theorem [2]) 

Since Vs(e) = e n ^ / u (e) and 

u&V 


/”( £ ) = 2riS(/ e 


/ r c + io ° adj(I — 4(s)) t 


det(/ — 4(s)) 


■M 


VgM 


v&V 


(s) £ s ds 











































Graph-Directed Sprays and Their Tube Volumes via Functional Equations 


22 


we have to evaluate the integral on the right hand side. As is well-known there is a 
general procedure to evaluate this integral by applying the residue theorem. In the 
present case of graph-directed sprays however, the det (I — A(s)) in the denominator 
calls for a more cautious treatment in order to be able to give an explicit validity range 
of the tube formula for small e. 

The above sum consists of integrals of the type 


1 

27ri 


pc+ioo 


adj (I - A(s)) uv g s 1 
det (I — A(s)) s — i 


e s ds (i = 0,1,..., n), 


for n — 1 < D < c < n. Recall that det (I — A(s)) is non-zero for Re(s) > D. 

We first note that there exists a vertical strip containing all the zeros of det (I—A(s)). 
To see this, notice that det (I — A(s)) can be expressed as a sum 


1 +J2p s o‘-J2^ 

a 0 

with 0 < p a ,qp < 1 since the entries of the matrix A(s) are of the form ^ r s e 
for 0 < r e < 1. As Re(s) —>• —oo, the smallest of p a ,qg will dominate and avoid 
det (/ — A(.s)) to vanish. We choose a q < 0 such that | det (I — A(s))| > 5 for some 
5 > 0 and for all Re(s) < q. 

We will choose a sequence —> oo such that | det(/ — A(s))| will be uniformly 

away from zero on the line segments ci < Re(s) < c, Im(s) = ±7j. 

Lemma 4 There exists an increasing sequence T j —t oo and a K > 0 such 

that | det (I — A(s))| > K for ci < Re(s) < c and Im(.s) = ±Tj. 

Proof. Being an entire function, det(/ — A(s)) has isolated zeros and we can choose 
7~i > 0 such that there are no zeros on the segment [q, c] x {n}. Let 2 K be the minimum 
of | det (I — A(s))| on this segment. 

To construct 72, we need the following lemma. 
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Lemma 5 (Dirichlet Lemma, |l]) Let M, N € N and T > 0. Let a i, cq,..., ajy be 
real numbers. There exists a real number h € [ T,T M N \ such that 


ai h\\ < 


1 

M 


(1 <i< N). 


Here || • || denotes “the distance to the nearest integer” function on R. 

The idea for choosing T2 = t\ + h will be the following. We want to arrange h such 
that 

| det(/ — A(s + ih)) — det(/ — A(s))| < K 

for q < Re(s) < c. Then the minimum of | det(/ — A(s + ih)) | will be greater than K. 
Since det (I — A(s)) = 1 + p s a — q^, we have 

a p 


det(/ — A(s + ih)) — det(/ — A(s))| = 


E 


< £ k + "'-Pc 


rC' h - K 


E 

P 


s+ih 


+ £ 


s-\-ih 




E* 


Re(s) 

a 


Pa-l 


+E«? 

s 


Re(s) 




Since the number of terms and Pq C ^ s \ q^ e ^ are bounded, it will be enough to make 
the factors |p^ 1 — 11 , q 1 ^ — 1 small enough. To realize this, we can apply the Dirichlet 
Lemma (Lemma 0 to make p 1 ^ = e lhlnpa and q l p = e lhlnq P close enough to 1, by 


choosing 


h In p a 


hlnqp 

2 vr 


2 vr 


small enough. 


As we have control on choosing h on any range, we can repeat this procedure to get 


a sequence of segments [q,c] x {rj} on all of which | det(/ — A(s))| is bounded below 
by K. m 


Let us now consider the rectangles Rj = [q,c] x [—Tj,Tj] and denote its oriented 
edges by L\ j. L 2 j, L 3 j, L^j as shown in Figure flOl We will show that for small enough 
e the integrals on L 2 j,L^j and L^j will tend to zero as j —» 00 so that by residue 
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theorem we will get the integral on the vertical line at c as a series of residues on the 
strip ci < Re(s) < c: 

1 [ c+i °° adj(I ~ A(s)) uv g s - 1 _ v /adj (I - A(s)) uv \ 

27ri./ c _ ioo det(J-A( a )) s-i £ ' ^ w re V det(I-A(s)) s-i’ U ) 

for i = 0,1 ,,n — 1. For i = n the same formula holds with the only difference that 
the residues are taken on 2 ). 



Figure 10: The strip containing the poles of the geometric zeta function C«( s ) and the 
rectangle Rj with the oriented boundary segments L i j , L 2 j, L 3 J, L^j used in the proof. 


First consider the integral on L 2 J = t + irj, ci < t < c. 


f adj(I — A(s)) uv g s 1 _ 

< 

r 

adj(/ - A(t + i Tj)) uv 

(/ 1 

Jl 2J det (I - A(s)) s-i £ 


J Cl 

det(7 — A(t + i.Tj)) 

\t + i Tj ■ 


adj(J — is of the form I + p s a — q s p (I indicates that 1 might be present 

a p 

or absent) and therefore is bounded for q < Re(s) < c, so that we can write 


f ad)(I - A(s)) uv g s 1 _ 

Jl 2J det(J-A(s)) s-i £ 


< 



(C being a constant not depending on j) 


which tends to zero for tj —>• 00 . 

Similarly the integral on L 4j - —>• 0 for Tj 00 . 
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Finally we consider the integral on L 3 j. 



adj(/ - g s 1 _ 

det(7 — .A(s)) s — i £ 



adj(J - A(s))^ g s 1 _ 
det (I — .A(s)) s — i £ 


where Cj = q + Tje 1 *, f < t < (see Figure fTTI) . adj (I — A(s)) uv is of the form 
1 + YlaPa — YlpQp- I n an y term, p a ,Q/3 there can appear at most (N — 1). power of 
the smallest weight r m i n of the graph where N is the number of the nodes of the graph. 
We can thus dominate |adj(J — A(s)) u „| by C We had | det(/ — v4(s))| > 5 

by the choice of q. 



Figure 11: The semi-circle Cj used to evaluate the integral on the segment L 3 


<3- 


We can now write 


&d](I - A(s)) uv g s 


ICi det (I - A(s)) s-i 


e s ds 


< 


L 


Q! r (N-l)Re(s) R e(s)-i 


Re(s) 


Cj u I s — l\ 

3tt / jy_^ ^ T~j COS t 


|ds| 


< C" 


< c" 


2 r ■ a 

I min & 


,37r / jy_^ Tj COS t 


| ci + Tje u — i\ 


dt 


2 r ■ a 

/ mm & 


dt 


since Tj < \ci + r^e 1 * — i\. By the Jordan Lemma (which states that 


» 3 zl 
2 


lim / a n cos dt = 0 , 

n—>-oo / 7r 
J 2 
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for any fixed a > 1) this integral tends to zero if 


r 


N -1 
min 


9 


£ 


> 1 . 


So we come to the conclusion that our tube formula is valid pointwise for 

e < r^T 1 min{g u } for monophase generators. In the pluriphase case, one should take 
uGV 

£ < r^ 1 min{gi )U }. Needless to say, these are only sufficient bounds. 
uGV 

Remark 7 There is an unfortunate misprint in Q p. 159, line 12], where £ < g 
should be £ < g\. 
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